Hub and spoke system for detecting and locating gas leaks

ABSTRACT

A system for detecting gas leaks and determining their location and size. A data gathering portion of the system utilizes a hub and spoke configuration to collect path-integrated spectroscopic data over multiple open paths around an area. A processing portion of the system applies a high resolution transport model together with meteorological data of the area to generate an influence function of possible leak locations on gas detector measurement paths, and applies an inversion model to the influence function and the spectroscopic data to generate gas source size and location.

This invention was made with government support under grant numberDE-AR0000539 awarded by the U.S. Department of Energy. The governmenthas certain rights in the invention.

U.S. Pat. No. 8,564,785 is incorporated herein by reference. Provisionalapplication for patent No. 62/160,163, filed 12 May 2015 is incorporatedherein by reference.

BACKGROUND OF THE INVENTION

Field of the Invention

The present invention relates to apparatus and methods for detecting andlocating gas leaks. In particular, the present invention relates to ahub and spoke spectroscopy system for detecting and locating methaneleaks.

Discussion of Related Art

Current open-path techniques that are capable of measuring methane leaksover long paths include diode laser-based absorption systems, LIDARsystems, and FTIR-based systems. Mobile FTIR systems suffer from lowwavelength resolution (large instrument distortion), and have thus faronly demonstrated ˜5-10% measurement uncertainty for trace GHGs, whichis far too great to detect, locate, and size methane leaks at kilometerscale standoff distances.

Diode laser-based systems and LIDAR systems focus on measurements of afew wavelengths around a single absorption feature of methane (or awavelength sweep over 1-2 features). High precision, long-termstability, and accuracy is difficult due to turbulence-induced laserintensity fluctuations and interference from overlapping absorption ofother molecules that are not included in spectral fits. Even techniqueswhich rely on detection of phase shifts induced by absorption features(instead of direct absorption) must account for phase shift induced byany absorbing component in the beam path and neighboring absorptionfeatures.

Sparse wavelength laser systems also do not typically measure otherspecies, temperature, pressure, or water vapor. A simultaneousmeasurement of water vapor, temperature, and pressure is desirable forcorrecting measured methane mole fractions to dry-air mole fractions, toaccount for time varying dilution effects of water vapor change on theapparent concentration of methane. In addition, water vapor, temperatureand pressure influence methane absorption feature shape, which isimportant when fitting the absorption features to accurately extract themethane mole fraction for calibration-free operation.

Many previous methane studies near oil and gas operations were performedwith commercial cavity-ringdown laser spectrometers (CRDS) either fixed,or mounted on vehicles and aircraft. These spectrometers enable veryhigh sensitivity with short measurement times, but require periodiccalibration, and are expensive. For specific leak detection withinversion techniques, the sensors either require an operator (pilot ordriver) or a network of multiple expensive sensors and commoncalibration.

Several other types of low-cost in-situ sensors for methane exist. Somefocus on making flux measurements because they are not stable over longperiods of time. Others lack the measurement precision needed toidentify smaller leaks or need to be calibrated often and corrected foreffects of temperature, pressure, humidity, or other interfering species(possibly requiring regular access to the well pad). Other in-situsensors with lower cost than CRDS sensors still require either anoperator to get spatial information or multiple sensors. Using multiplesensors requires inter-calibration and inter-comparability between thevarious sensors to correct for background fluctuations in methane with aremote background sensor or to compare methane concentration betweensites. In a distributed system, each sensor may require power andcommunication.

A need remains in the art for apparatus and methods for detecting gasleaks capable of sensitivity, accuracy, lack of calibration, andmulti-species operation over kilometer-scale paths.

SUMMARY OF THE INVENTION

It is an object of the present invention to provide apparatus andmethods for detecting gas leaks capable of sensitivity, accuracy,calibration-free operation, and multi-species detection overkilometer-scale paths.

A low cost dual comb spectrometer design uses the method of deploying aline of sight, broadband, laser absorption sensor to locate and sizetrace gas leaks. An embodiment of the present invention includes an openpath spectrometer gas detector combined with a high resolution transportmodel and an inversion model to detect and locate gas leaks.

A system according to the present invention determines the location andsize of a gas source within an area by providing a spectrometer gasdetector, collecting path-integrated spectroscopic data over multipleopen paths around the area with the detector, collecting meteorologicaldata related to the area, applying a high resolution transport modeltogether with meteorological data to generate an influence function ofpotential source locations on gas detector measurement paths, andapplying an inversion model to the influence function and spectroscopicdata to generate gas source size and location. The spectrometer might bebased on a dual comb spectrometer. In a preferred embodiment, the highresolution transport model is a large eddy simulation. The system mightemploy a number of reflectors, such as retroreflectors arrayed aroundthe area, or a mobile reflector, for example on a UAV. Themeteorological data might be measured onsite or provided by a simulationmodel of the area. Some embodiments use Kalman filtering to update theinversion model.

The system might include a telescope for transmitting laser beams andreceiving the reflected beams, and a gimbal for orienting the telescopeto scan the area.

BRIEF DESCRIPTION OF THE DRAWINGS

FIGS. 1A-1C (Prior Art) illustrate the operation of a dual combspectroscopy (DCS) system that is useful in the data gathering portionof the present invention. 1A (Prior Art) is a schematic block diagramillustrating the process of passing two frequency combs through a gasand detecting the resulting light. FIG. 1B (Prior Art) shows the twofrequency combs after the light passes through the gas. FIG. 1C (PriorArt) illustrates the resulting heterodyne interference signals.

FIG. 2 is a diagram illustrating an embodiment of the data gatheringportion of the present invention.

FIG. 3 is a block diagram illustrating an embodiment of the entiresystem according to the present invention.

FIG. 4A is a flow diagram illustrating a general embodiment of thesignal processing portion of the present invention. FIG. 4B is a flowdiagram illustrating a specific embodiment of the preferred signalprocessing portion of the present invention.

DETAILED DESCRIPTION OF THE INVENTION

FIGS. 1A-1C (Prior Art) illustrate the operation of a dual combspectroscopy (DCS) system that is useful in the gas detection portion ofthe present invention 1A (Prior Art) is a schematic block diagramillustrating the process of passing two frequency combs (e.g. nearinfrared light) having slightly different tooth spacings through a gasand detecting the resulting light. FIG. 1B shows the two frequency combsafter the light passes through the gas, so that some light frequencieshave been absorbed by the gas. FIG. 1C illustrates the resultingheterodyne interference signals at the detector, for example at rffrequencies.

DCS overcomes the key limitations of single or sparse wavelengthabsorption or LIDAR approaches: it enables accurate correction of thebaseline laser intensity, and simultaneous measurement of CH₄, ¹³CH₄,H₂O, other species (such as ethane and propane), temperature, andpressure. It therefore reports interference-free, true dry-air molefractions that account for variable water vapor dilution. With noinstrument distortion (of lineshape) and a near perfect wavelength axis,the technique is also drift-free and requires no calibration. Comparedwith single point measurements that might be deployed on a tower ormobile platforms (aircraft or cars), this solution requires no operatorinvolvement and can interrogate multiple locations simultaneously.

FIG. 2 is a diagram illustrating an embodiment of the data gatheringportion 200 of the present invention. The open path hub-and-spoke systemincludes a central spectrometer/detector unit 202 and several reflectors210 arrayed over an area to allow unit 202 to transmit and receive lightin a variety of directions to detect gas leaks 208 from a number ofwells 206. Unit 202 scans the area with laser beams 204 over long openpaths and detects the light reflected from reflectors 210. Data isgenerally collected over a period of days under various weatherconditions.

A processor 330 (see FIGS. 3 and 4) combines the data from the reflectedbeams 328 with current meteorological data 340 (such as the directionand speed of wind 220) of the area to determine the location of any gasleaks. For example, processor 330 combines data from reflected beamswith high-resolution computational fluid dynamics and inversiontechniques to locate and size the leaks.

In some embodiments, reflectors 210 are retroreflectors which reflectthe light directly back to the spectrometer/detector unit 202.Retroreflectors provide a high degree of pointing flexibility (the beamsfrom a large array of incoming angles are re-directed back the directionthey came in), and thus alignment is automatic as long as the beam cantrack the retro-reflector.

Reflectors 210 might be located on well platforms or risers, towers,trees, fences, etc. Alternative reflectors may be used, including otherfixed reflectors and even environmental reflectors such as buildings. Insome embodiments, a UAV 212 with a reflector 214 periodically flies afixed path around the area and beams 204 are reflected off of reflector214 in various directions covering the area containing wells 206.

One particularly useful embodiment includes a dual comb spectrometer(DCS) unit 300 (See FIG. 3). Since DCS 300 is compatible with fiberoptic transmission, it is possible to monitor the perimeter of largefacilities using several paths serviced by a single spectrometer viafiber optics. The same is true of multiple centralized towers separatedby many kilometers. A single DCS spectrometer can supply light viaoptical fiber to multiple towers to reduce system cost.

In a preferred embodiment, a central, high-performancespectrometer/detector unit 202 sends light 204 sequentially overdifferent long open paths to retroreflectors 210. The measuredabsorption spectra are fit to determine methane concentrationenhancements. High resolution Large Eddy Simulation-based inversiontechniques are used to interpret the measured enhancement into methaneleak size and location. The 1 ppb sensitivity of spectrometer/detectorunit 202 enables detection of small leaks over a range of heights anddownwind distances, thus enabling flexible location of retroreflectors210. Accurate, sensitive, calibration- and drift-free measurements ofmethane have been demonstrated over a 2 km open air path.

FIG. 3 is a block diagram illustrating an embodiment of the entiresystem according to the present invention. Spectrometer/detector unit202 transmits DCS beams 204 via transmitter/receiver 322 towardreflectors 210 and detects the reflected beams. Processor 330 combinesdata 328 from the reflected beams with meteorological data 340 to detectand locate any gas leaks (see FIG. 4). Meteorological data is needed tomodel plume shape and path of gas leaks 208.

DCS unit 300 includes two comb units 302 and 304, as well as circuitry306 to stabilize the generated comb combination 308 and/or toelectronically post-correct the spectroscopic signal 328. Combs 308 areprovided to transmitter/receiver unit 322 including telescope 326 fortransmitting and receiving light beams 104, gimbal 324 to scan beams 204over the area, and detector 320 to detect the reflected beams fromreflectors 210 and provide data 328 to processor 330.

Processor 330 provides control signals 332, 334 to DCS unit 300 andtransmitter/receiver unit 322. It also receives reflected beam data 328from detector 320 and meteorological data 340. Processor 330 outputs 350the location of any gas leaks within the area.

Meteorological data 340 might comprise local measurements of wind 220,as well as temperature, humidity, etc. Alternatively, it could bedetermined from meteorological simulations of the area, such as WRF. Ora combination of these methods may be used.

FIG. 4A is a flow diagram illustrating a general embodiment of thesignal processing portion 330 of the present invention. Meteorologicaldata 340 are provided to a high resolution transport model 402. Highresolution transport model 402 is used to create an influence functionthat is passed together with spectrometer data 328 to the inversionmodel 404 that computes the leak location and size 350.

Standard kilometer-scale atmospheric transport models such as thewidely-used Weather Research Forecast (WRF) model may not have desiredresolution to locate a target leak on the meter-scale that is requiredby the MONITOR program. Spectrometer data can be coupled withcomputational fluid dynamics (CFD) model of the atmospheric transportand an inversion technique to locate and size methane leaks.

FIG. 4B is a flow diagram illustrating a specific embodiment of thepreferred signal processing portion of the present invention. Highresolution large eddy simulations (LES) 402 incorporate the localtopography (elevation, large buildings, vegetation, etc.) and localmeasurements of wind, temperature, and pressure for the inversiontechnique. Meteorological data 340 is provided as an input to the LargeEddy Simulations 402. An adjoint of the LES model is computed todetermine sensitivities. These sensitivities are combined withspectrometer data as well as prior leak information to obtain updatedleak information and uncertainties as part of a Kalman filteringapproach 404.

The NIST Fire Dynamics Simulator (FDS) is a high-resolution CFD toolthat was recently demonstrated for transport of CO₂ in an urbanenvironment. Originally developed for fire-driven flows, NIST's FireDynamics Simulator (FDS Version 5.5) is an open-source CFD codecontaining a hydrodynamic solver suitable for low-speed flow (Mach <0.3)and a thermal radiation transport model.

FDS is used mainly for large eddy simulations (LES), in which fluidmotion and turbulent eddies are resolved at scales greater than the meshcell size, while the dissipative effects of turbulence at sub-gridscales are modeled with a turbulent viscosity, rather than directnumerical simulation, which would require ˜10¹⁵ cells in a simulationdomain only 100 m on a side. The LES technique might employ cells of 5 maverage dimension, which requires about 7×10⁷ cells to model theatmospheric boundary layer in a domain with 3 km radius.

While the exemplary preferred embodiments of the present invention aredescribed herein with particularity, those skilled in the art willappreciate various changes, additions, and applications other than thosespecifically mentioned, which are within the spirit of this invention.For example, this technique could be used to locate and determine theflux of any gas source. The technique could be extended to look forother trace gasses besides methane to, for example, look for chemicalleaks in industrial facilities or to, for example, look for CO₂ leaks ata carbon sequestration site. Other extensions could determine emissionsfrom industrial sites, agricultural sites, animal-raising operations, orchemical and biological weapons releases.

What is claimed is:
 1. A method of determining the location and size ofa gas source within an area comprising the steps of: collecting, using aspectrometer gas detector, path-integrated spectroscopic data overmultiple open paths around the area; applying a high resolutiontransport model together with meteorological data to generate aninfluence function of potential source locations on gas detectormeasurement paths; determining gas source size and location by applyingan inversion model to the influence function and spectroscopic data. 2.The method of claim 1, the spectrometer gas detector being a dual combspectrometer.
 3. The method of claim 2 the high resolution transportmodel being a large eddy simulation.
 4. The method of claim 1 the highresolution transport model being a large eddy simulation.
 5. The methodof claim 1, the step of collecting path integrated spectroscopic dataincluding reflecting light from the spectrometer gas detector off ofmultiple retroreflectors and detects the reflected light with areceiver.
 6. The method of claim 1, the step of collecting pathintegrated spectroscopic data including reflecting light from thespectrometer gas detector off of a mobile reflector at differentlocations and detecting the reflected light with a receiver.
 7. Themethod of claim 1, the meteorological data being measured from withinthe area.
 8. The method of claim 1, the meteorological data being from ameteorological simulation of the area.
 9. The method of claim 1 furtherincluding the step of using Kalman filtering to update the inversionmodel.
 10. Apparatus for determining the location and size of a gassource within an area comprising: a spectrometer unit having a beamtransmitter and a detector, the unit configured to scan the area andcollect path-integrated spectroscopic data; multiple reflection pointswithin the area for reflecting beams from the transmitter such that theyare detected by the detector; a source of meteorological data for thearea; a processor configured to receive data based on the detected beamsand the meteorological data, to apply a high resolution transport modelto the meteorological data; and to apply an inversion model to the highresolution transport model and the spectroscopic data; therebycalculating location and size of a gas source.
 11. The apparatus ofclaim 10 wherein the spectrometer unit includes a dual combspectrometer.
 12. The apparatus of claim 10 wherein the multiplereflection points comprise multiple reflectors arrayed around the area.13. The apparatus of claim 10 wherein the multiple reflection pointscomprise a mobile reflector configured to move around the area.
 14. Theapparatus of claim 10 wherein the spectrometer unit includes a telescopefor transmitting laser beams and receiving the reflected beams.
 15. Theapparatus of claim 14, further comprising a gimbal for orienting thetelescope.
 16. Apparatus for determining the location and size of a gassource within an area comprising: means for collecting path-integratedspectroscopic data over multiple open paths around the area; means forcollecting meteorological data related to the area; means for applying ahigh resolution transport model together with meteorological data togenerate an influence function of potential source locations on the openpaths; means for applying an inversion model to the influence functionand spectroscopic data to generate gas source size and location.
 17. Theapparatus of claim 16 wherein the means for collecting includes a dualcomb spectrometer.
 18. The method of claim 17 wherein the highresolution transport model comprises a large eddy simulation.
 19. Themethod of claim 16 wherein the high resolution transport model comprisesa large eddy simulation.
 20. The method of claim 16 further includingKalman filtering means for updating the inversion model.